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A method is proposed to characterize and quantify multipartite entanglement in terms 
of the probability density function of bipartite entanglement over all possible balanced 
bipartitions of an ensemble of qubits. The method is tested on a class of random pure 
states. 

The quantification of multipartite entanglement is an open and very challenging 
problem. An exhaustive definition of bipartite entanglement exists and hinges upon 
the von Neumann entropy and the entanglement of formation, 1 but the problem of 
defining multipartite entanglement is more difficult 2 and no unique definition exists: 
different definitions tend indeed to focus on different aspects of the problem, cap- 
turing different features of entanglement, 3 that do not necessarily agree with each 
other. Moreover, as the size of the system increases, the number of measures (i.e. 
real numbers) needed to quantify multipartite entanglement grows exponentially. 

This work is motivated by the idea that a good definition of multipartite entan- 
glement should stem from some statistical information about the system. 4 We shall 
therefore look at the distribution of the purity of a subsystem over all bipartitions 
of the total system. As a measure of multipartite entanglement we will take a whole 
function: the probability density of bipartite entanglement between any two parts 
of the total system. According to our definition multipartite entanglement is large 
when bipartite entanglement (i) is large and (ii) does not depend on the bipartition, 
namely when (i+ii) the probability density of bipartite entanglement is a narrow 
function centered at a large value. This definition will be tested on two class of 
states that are known to be characterized by a large entanglement. We emphasize 
that the idea that complicated phenomena cannot be "summarized" in a single (or 
a few) number(s) was already proposed in the context of complex systems 5 and has 
been also considered in relation to quantum entanglement. 6 

We shall focus on a collection of n qubits and consider a partition in two sub- 
systems A and B, made up of tia and Ub qubits (tia + tlb = n), respectively. For 
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definiteness we assume ha < n B . The total Hilbert space is the tensor product 
H = H A ® H B and the dimensions are dimW = A = 2™, dimH A = Na = 2" A and 
dimW B = 7V B = 2" B , respectively (N A N B = A). 

We shall consider pure states. Their expression adapted to the bipartition reads 

N-l N A -1N B -1 

W) = E zfc i fc > = E E ^jbIja) ® |/ B >, (i) 

fe=o i^=o i B =o 

where \k) = Ij'a)® Kb), with a bijection between fc and (ja,Ib)- Think of the binary 
expressions of an integer k in terms of the binary expression of {]a,Ib)- 

As a measure of bipartite entanglement between A and B we consider the par- 
ticipation number 

Nab = n A B> ttab(I^)) = tr^pi, Pa = tr B p, (2) 

where p = \ip)(ip\ and tr^ (trg) is the partial trace over the subsystem A (B). 
Nab measures the effective rank of the matrix pa, namely the effective Schmidt 
number. 7 We note that 

1 < N AB = N BA < min(JV A , N B ), (3) 

with the maximum (minimum) value attained for a completely mixed (pure) state 
Pa- Therefore, a larger value of Nab corresponds to a more entangled bipartition 
(A,B), the maximum value being attainable for a balanced bipartition, i.e. when 
ha = [n/2] (and n B = [(n + l)/2]), where [x] is the integer part of the real x, 
that is the largest integer not exceeding x. The maximum possible entanglement is 
Nab = Na = 2™ A . The quantity uab = log 2 Nab represents the effective number 
of entangled qubits, given the bipartition (namely, the number of bipartite entan- 
glement "links" that are "severed" when the system is bipartitioned). 

Clearly, the quantity Nab will depend on the bipartition, as in general entan- 
glement will be distributed in a different way among all possible bipartitions. As 
explained in the introduction, we are motivated by the idea that the distribution 
p(Nab) of Nab yields information about multipartite entanglement. 

Let us therefore study the typical form of our measure of multipartite entan- 
glement p(Nab) for a very large class of pure states, sampled according to a given 
symmetric distribution on the projective Hilbert space {ip € Tt, = 1} (e.g. the 
unitarily invariant Haar measure). By plugging (1) into (2) one gets 



kab 



N A -\ N 

E 



E Z 3l Z 3'l Z 3'l' Z jl' ( 4 ) 

j,j'=o U'=o 

and it can be shown 4 that, in the thcrmodynamical limit (that is practically attained 
for n > 5), independently of the distribution of the coefficients, the mean and the 
standard deviation of (4) over all possible balanced bipartitions read 

Na + N b -1 fa 2 2 
Pab = ^ = V A ' ° AB = N*' ( ge > ( > 
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respectively, where a = 8/2 (a = 9/2) for even (odd) n. Moreover, the probability 
density of Nab m Eq. (2) reads 

It is interesting to compare the features of these generic random states with those of 
other states studied in the literature. Table 1 displays the average value of Nab for 

Table 1. Mean bipartite entanglement (Nab) for 
different states and different number of qubits n. 



n 


GHZ 


W 


cluster 


random 


5 


2 


1.923 


3.6 


2.909 


6 


2 


2 


5.4 


4.267 


7 


2 


1.96 


6.171 


5.565 


8 


2 


2 


8.743 


8.258 


9 


2 


1.976 


10.349 


10.894 


10 


2 


2 


14.206 


16.254 


11 


2 


1.984 


17.176 


21.558 


12 


2 


2 


23.156 


32.252 



GHZ states, 8 W states, 9 cluster states 10 and the generic states (1), for n = 5 -j- 12. 
While the entanglement of the GHZ and W states is essentially independent of n, 
the situation is drastically different for cluster and random states. In both cases, 
the average entanglement increases with n; for n > 8 the average entanglement is 
higher for random states. However, the mean (Nab) yields poor information on 
multipartite entanglement. For this reason, it is useful to analyze the distribution 
of bipartite entanglement over all possible balanced bipartitions. 

The results for the cluster and random states are shown in Fig. 1, for n rang- 
ing between 5 and 12. Notice that the distribution function of the random state 
is always peaked around (Nab) — Hab given by (5). Notice also that the cluster 
state can reach higher values of Nab (the maximum possible value being 2^ n / 2 ^), 
however, the fraction of bipartitions yielding this result becomes smaller for higher 
n. This is immediately understood if one realizes that the cluster states are designed 
for optimized applications and therefore perform better in terms of specific biparti- 
tions. On the other hand, according to the measure we propose, the random states 
are characterized by a large value of multipartite entanglement, that is roughly 
independent of the bipartition. 

In Fig. 2 we compare the number of balanced bipartitions vs Nab for the random 
states and increasing n. The related probability density functions (6) are displayed 
in Fig. 3. Notice that as the number of spins increases from n — 5 to n — 12 the 
mean increases and the distribution becomes relatively narrower. As we emphasized, 
these are both signatures of a very high degree of multipartite entanglement, whose 
features become (as n increases) practically independent of the bipartition. In Fig. 
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P n p p n p 




2 3 4 5 6 7 8 Nab 2.5 5 7.5 10 12.5 15 N A B 

pn p p n p 




5 10 15 20 25 30 Nab 10 20 30 40 50 60 Nab 

Fig. 1. Number of balanced bipartitions vs Nj^g\ n v = n\/n^\ng\ is the number of bipartitions. 
the light-gray bars represent cluster states, the dark-gray ones random states; the solid line is the 
distribution (5)-(6); the black arrows indicate the average (N^g) c i uateI . For even n (n = 12 in 
particular) the distribution of the random state partially hides a peak of the corresponding cluster 
state distribution, centered at N AB = 2 n A~ 1 = 2I™/ 2 !- 1 . 
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3 it is interesting to observe the difference between the distributions for odd and 
even n. 




2 4 6 8 10 Nab 10 15 20 25 30 N AB 

Fig. 2. Number of balanced bipartitions vs Nab\ the histograms are numerically obtained for 
the typical states; the solid line represents their distribution. The value of n (total number of 
spins) is always indicated and ranges (a) from 5 to 8; (b) from 9 to 12; n v = nl/nAlng] is the 
number of bipartitions. 
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Fig. 3. Probability densities functions (6) vs Nab ■ Each curve is labeled with the corresponding 
value of n (number of qubits). The standard deviation c quickly becomes independent of n [sec 
Fig. 4(b)] and depends only on parity of the latter. 



In Fig. 4(a) we plot the value of (Nab) for the cluster and random states (see 
Table 1). We notice that, for n = 9, (Nab) random becomes larger than (iV J 4s) c i ust or- 
Figure 4(b) displays the behavior of the standard deviation of Nab, 

ct~(j A b/ii 2 ab . (7) 

For the cluster states this quantity tends to diverge when the size of the system in- 
creases. By contrast, from Eq. (5), o = \/2/a is constant for the typical states. This 
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Fig. 4. Comparison between typical and cluster states: (a) expectation value {Nab) and (b) 
standard deviation a of the distributions in Fig. 3 vs n (number of qubits). 



means that the ratio a /{Nab} tends to 0. Finally, Figure 5 displays a parametric 
plot of a vs (Nab)- Clearly, for the random states a is independent of (Nab)- 
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* * * * * * * * 

5 10 15 20 25 30 (N AB ) 



Fig. 5. Standard deviation <r vs expectation value (Nab) of the distributions in Fig. 3. 



We emphasize that our analysis should by no means be taken as an argument 
against the performance of the cluster states. As we stressed before, cluster states 
are tailored for specific purposes in quantum information processing, and in that 
respect are very well suited. We compared the generic states to the cluster states 
specifically because the latter are also known to be characterized by a large entan- 
glement. 

An efficient way to generate states endowed with random features is by means 
of a chaotic dynamics, 12 or at the onset of a quantum phase transition. 13 In partic- 
ular, the random states describe quite well states with support on chaotic regions 
of phase space, before dynamical localization has taken place. These features make 
these states rather appealing, from a practical point of view, in that they are eas- 
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ily generated. The introduction of a probability density function as a measure of 
multipartite entanglement paves the way to further investigations of the intimate 
relation between entanglement and randomness and their behavior across a phase 
transition. 
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